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REMOTELY  SENSED  WAVE  SPECTRA  FROM 
JOINT  NORTH  SEA  EXPERIMENTS  II 

I . INTRODUCTION 

The  Joint  North  Sea  Wave  Project  II  (JONSWAP  II)  was  a 
cooperative  effort  by  a number  of  scientists  from  Germany, 
Holland,  Denmark,  England,  Canada  and  the  United  States. 

One  of  the  experimental  objectives  was  to  obtain  remote 
sensing  data,  from  an  instrumented  aircraft  simultaneously 
with  the  ground-truth  data  from  an  array  of  surface  instru- 
ments for  addressing  important  questions  regarding  the 
interpretation  of  remotely  sensed  data. 

Scientists  from  the  Naval  Research  Laboratory  and  NASA 
Wallops  Flight  Center  utilized  NASA's  instrumented  C-54 
aircraft  for  observation  in  JONSWAP  II,  as  one  of  the  remote 
sensing  teams.  A typical  remote  sensing  instrument  carried 
aboard  the  aircraft  was  the  laser  profilometer.  This  air- 
borne laser  profilometer  profiled  the  ocean  surface  as  the 
aircraft  flew  with  constant  velocity  and  at  the  altitude  of 
91.5  or  152.4  meters  above  the  ocean.  Wave  spectra  calcu- 
lated from  wave  data  obtained  on  September  19,  1973,  are 
presented  in  this  report.  The  data  were  relatively  less 
noisy  and  steady  offshore  wind  conditions  were  observed  on 
that  particular  day. 

Due  to  the  lack  of  strong  current  gradients  in  the 
experiment  area,  the  original  goal  of  demonstrating  the 
feasibility  of  inferring  currents  by  using  the  technique 
proposed  by  Huang,  et  al.  (1972)  could  not  be  properly 

Note:  Manuscript  submitted  June  15,  1978. 
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investigated.  Nevertheless,  the  technique  has  been  worked 
out  and  included  in  this  report  as  an  Appendix. 


The  placements  of  in-situ  ground-truth  instruments  of 
JONSWAP  II  are  presented  in  Figure  1.  The  instruments  at 
Stations  4 and  8 are  shown,  particularly,  for  their  arrange 
ments.  Airborne  measurements  were  made  in  the  general  area 
from  Station  4 to  Station  11  during  the  period  of  September 
17  to  September  27.  Of  which,  only  September  18  and 
September  19  had  the  condition  of  steady  offshore  wind  in 
which  the  analysis  is  interested  because  of  its  assumptions 
the  data  taken  on  these  two  days  is,  thus,  selected  for  the 
study. 


The  flight  patterns  of  these  two  days  are  shown  in 
Figures  2 and  3.  According  to  Dr.  Walsh  (private  communi- 
cation), of  NASA  Wallops  Flight  Center,  the  NASA  airplane 
was  flying  at  the  altitude  of  152.4  meters  on  September  18 
and  at  the  altitude  of  91.5  meters  on  September  19.  The 
small  circles  on  the  flight  paths  in  Figures  2 and  3 indicate 
the  locations  at  which  data  from  Inertial  Navigation  System 
was  taken.  The  symbol  "P"  indicates  that  photographs  were 
taken  with  70  mm  Hasselbold  cameras  (80  mm  lens)  out  the 
left  side  of  the  aircraft  at  the  depression  angles  of  45 
degrees . 


On  September  18  over  Station  8 the  wind  was  recorded  as 
steady  for  14  hours  at  13.5  M/sec  in  the  direction  of  121 
deg.  from  the  north.  While  on  September  19  over  the  same 
station  the  wind  was  in  the  direction  of  147  deg.  from  the 
north  and  remained  steady  for  12  hours  at  12.5  M/sec.  The 
detailed  information  on  these  two  flights  as  well  as  the 
available  in-situ  measurements  is  listed  in  Table  1.  The 
dotted  lines  indicate  that  the  information  was  not  recorded 


Table  1 - Airborne  and  In-Sltu  instruments 


but  can  be  estimated.  The  heading  of  the  airplane  is 
presented  in  a general  sense  of  upwind  and  downwind,  however, 
this  information  is  used  precisely  in  the  calculation  as 
what  has  been  recorded.  Three  tapes  were  recorded  on 
September  18  and  four  tapes  were  recorded  the  next  day. 
Excluding  calibrations,  the  files  with  data  of  real  measure- 
ments are  listed  also  in  the  table.  Each  file  takes  91 
seconds  of  digitized  data  at  a frequency  of  90  Hz.  The  in- 
situ  instruments  listed  in  the  table  provide  wave  informa- 
tion. The  Institute  of  Oceanographic  Sciences  at  Wormley, 
Great  Britain,  operated  the  Telemetering  Accelerometer  Buoys 
A^  and  A2  at  Stations  10W  and  10E,  respectively.  The  Wave 
Rider  Buoys,  W^,  W^,  and  Wg  were  operated  by  the  Deutsches 
Hydrographic  Institute  at  Hamburg,  Federal  Republic  of 
Germany.  Very  unfortunately,  due  to  circumstances  beyond 
control  at  the  time  of  the  flights,  all  of  these  in-situ 
instruments  provide  a very  limited  quantity  of  data.  As  a 
consequence,  no  ground-truth  information  on  waves  at  Stations 
10  and  11  is  available  for  the  purpose  of  comparison. 

III.  RESULTS 

The  Spectra-Physics  Geodolite  3A  Laser  Profilometer 
used  in  this  experiment  happened  to  produce  very  noisy  data 
for  these  two  flights.  There  are  also  clear  indications 
from  the  data  that  this  laser  profilometer  lost  its  lock  on 
the  phase  shifting  mechanism  quite  frequently.  Among  all 
the  available  data  for  these  two  days,  data  from  File  4 of 
Tape  1 and  Files  1,  2,  and  3 of  Tape  2,  taken  on  September 
19,  is  barely  acceptable  in  quality.  After  careful  examina- 
tion these  four  files  were  chosen  for  further  analyses. 

Each  file  of  data  contains  8192  digitized  numerical 
values  at  the  time  interval  of  one-ninetieth  of  a second. 

A high  pass  filter  was  used  to  filter  out  the  effects  of 
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aircraft  motions  on  the  measurements.  The  numerical  filter 
applied  was  the  so-called  "Martin  Filter"  (Martin,  1957). 

It  is  a symmetric,  non-recursive  filter  which  features  a 
sine  termination  to  the  gain  function  and  a correction  which 
insures  unity  at  a cutoff  frequency  of  zero.  The  sine 
termination  is  introduced  to  avoid  large  oscillations  in  the 
gain  function  due  to  a sharp  cutoff.  The  cutoff  frequency 
used  in  the  filter  corresponded  to  a cutoff  period  of 
10  seconds  in  this  particular  application. 

Due  to  the  limited  storage  capacity  of  the  CDC  1604 
computer  this  amount  of  data  was  segmented  at  a time  length 
of  5.69  seconds  each.  This  length  in  time  is  considered 
adequate  for  obtaining  the  ocean  wave  spectrum.  Fifteen 
spectra  were  then  ensemble  averaged.  The  effects  due  to 
Doppler  shift  were  removed  by  assuming  that  waves  were 
travelling  in  the  direction  of  the  wind. 

Figures  4,  5,  6,  and  7 are  the  wave  spectra  calculated 

from  File  4 of  Tape  1,  and  Files  1,  2,  and  3 of  Tape  2, 

respectively.  In  addition.  Figures  8,  9,  10,  and  11  are  the 

wave  height  distribution  functions  obtained  from  the  files 

in  the  same  sequence.  The  vertical  coordinates  of  Figures 

4,  5,  6 and  7 are  in  logarithmic  scale  and  in  the  unit  of 
2 

M /Hz.  The  horizontal  coordinates  of  these  figures  are  also 
in  logarithmic  scale,  but  in  the  unit  of  Hz.  Figures  4 
through  7 show  wind  speeds  in  M/sec  and  fetches  in  M as  well 
as  Figures  8 through  11.  The  significant  wave  heights 
(SWH's)  in  M as  defined  by  Neumann  and  Pierson  (1966)  are 
alsq  shown  in  Figures  4 through  7.  The  abscissas,  of  Figure 
8 to  11,  represent  wave  amplitudes  and  are  nondimensionalized 
by  being  divided  with  variances  of  the  corresponding  wave 
spectra. 

The  laser  noise  levels  vary  from  0.048  M to  0.094  M 
while  the  laser  noise  limits  are  greater  than  3.47,  4.26, 


FREQUENCY  Hz 
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spectra  from  file  1,  tape  2,  of  September  19,  1973 
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MULTIPLES  OF  .784  METERS 
Wave  Amplitude 

Wave  height  distributions  from  file  4,  tape  1 
of  September  19,  1973 
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MULTIPLES  OF  .840  METERS 
Wave  Amplitude 

9 - Wave  height  distributions  from  file  1 
of  September  19,  1973 


3.36,  and  4.02  Hz  shown  by  Figures  4,  5,  6,  and  7,  respec- 
tively. These  are  unusually  high  noise  levels  at  low 
ranges  of  frequencies  which  are  the  indications  of  highly 
noisy  and  undesireable  laser  signals.  The  spectral  estimates 
which  appear  to  be  lower  than  laser  noise  levels  in  values 
at  the  low  frequency  are  caused  by  the  high  pass  filter. 

The  slopes  at  the  high  wave  frequency  ends  of  the  wave 
spectra  show  a range  of  values  from  -3.25  to  -4.71  in  Figures 
4 to  7.  Some  mechanisms  other  than  the  wave-breaking 
mechanism  suggested  by  Phillips  (1958)  and  Kitaigorodskii 
(1961)  must  have  been  in  control  of  the  physical  processes 
as  well  at  these  high  wave  frequency  ends.  For  otherwise 
the  slope  at  the  high  wave  frequency  end  should  show  a value 
of  5 if  the  wave-breaking  mechanism  were  the  only  one  which 
controls  the  growth  of  these  waves.  However,  due  to  the 
lack  of  strong  current  gradients  in  the  area  as  indicated  by 
the  in-situ  measurements  there  might  be  some  mechanisms 
other  than  wave-current  interaction  which  were  responsible 
for  the  variations  of  slope  values  at  the  high  wave  fre- 
quency ends  from  the  expected  value  of  “5. 

Nevertheless,  the  wave  height  distributions  shown  in 
Figures  8 to  11  display  the  characteristics  of  Gaussian 
distribution.  The  values  shown  beyond  2.1  M,  2.03  M,  and 
1.95  M in  Figures  9,  10,  and  11,  respectively  are  due  to  the 
drops  of  bits  in  the  recording  tapes  which  are  not  the 
results  of  the  real  physical  processes  in  the  ocean. 

The  remotely  sensed  significant  wave  heights  (SWH's), 
which  are  also  shown  in  Figures  4 to  7 and  are  defined  as 
four  times  the  variances  of  the  wave  spectra  (Neumann  and 
Pierson,  1966) , agree  well  with  those  from  visual  estimates 
(Walsh,  private  communications). 


IV.  CONCLUSION 


As  stated  by  Hasselmann,  et  al.  (1973)  for  waves 
generated  by  offshore  wind  at  the  site  of  JONSWAP  II,  the 
assumption  of  deep  water  is  valid  only  for  wind  speeds  below 
25  M/sec  and  the  assumption  of  homogeneous  wave  field  in  the 
direction  parallel  to  the  shoreline  is  valid  only  when  the 
angle  between  the  shoreline  and  wind  direction  is  less  than 
30  degrees.  Both  these  assumptions  are  valid  for  the  wind 
conditions  on  September  19,  upon  which  the  theoretical 
analyses  were  based. 

However,  due  to  the  lack  of  strong  current  gradients  in 
the  area  where  the  experiment  was  performed,  no  appropriate 
investigation  can  be  done  on  the  feasibility  of  inferring 
currents  by  employing  the  technique  proposed  by  Huang,  et 
al.  (1972).  Nevertheless,  this  proposed  technique  has  been 
worked  out  in  detail  and  is  included  in  this  report  as 
Appendix  A.  An  area,  with  known  strong  shear  current,  such 
as  the  western  boundary  of  the  Gulf  Stream  is  desireable  to 
be  established  as  the  test  site  for  inferring  currents  by 
employing  this  particular  technique. 

The  remotely  sensed  significant  wave  heights  (SWHs) 
agree  well  with  what  would  be  expected  for  SWHs  at  those 
particular  areas  in  the  North  Sea  at  the  time  of  the  experi- 
ments. 
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APPENDIX  A 


NUMERICAL  EVALUATION  OF  THE  INTERACTIONS 
BETWEEN  STEADY  NON-UNIFORM  CURRENTS  AND 
FETCH-LIMITED  WIND-GENERATED  GRAVITY  WAVES 


Al . THEORY 


a.  Theoretical  Background 

Let  n be  the  total  wave  frequency,  g be  the  gravita- 
tional acceleration,  U be  the  current  speed,  and  the  function 

R (n)  = (1  + 4Un/g)1/2  . (1) 

Huang,  et  al.  (1972)  have  shown  that,  under  the  influence  of 
a steady  non-uniform  current,  the  wave  frequency  spectrum, 
4>(n),  for  random  gravity  waves  in  deep  water  and  under 
steady  wind  conditions  is  given  by 

44>  (n) 

4>(n)  = t 

[1  + R (n)  ] [R  (n)  + (R(n)ri 

where  4>o  (n)  is  the  wave  frequency  spectrum  without  the 

influence  of  the  current.  Because  of  the  special  locality 

of  JONSWAP  II,  4>  (n)  is  chosen  to  be  the  fetch-limited  wave 

o 

frequency  spectrum  for  wind-generated  waves  as  proposed  by 
Hasselmann,  et  al.  (1973).  More  detailed  information  on 
this  wave  spectrum  will  be  presented  later  in  this  section. 

From  Equation  (2) , the  standard  deviation,  a,  of  wave 
amplitude,  n,  measured  from  mean  sea  level,  can  be  evaluated 
as 

(3) 

frequency 
30  cm. 
steady 
known . 


a2  = 


4)  (n)dn 


in  which  n is  the  cutoff  wave  frequency  of  the  wave 
c 

spectrum  taken  as  that  of  a wave  with  wave  length  of 
Thus,  from  Equations  (2)  and  (3),  the  influence  of  a 
non-uniform  current  on  the  standard  deviation,  a,  is 


A-l 


Subsequently,  as  done  by  Huang  (1974),  it  is  possible  to 
evaluate  the  probability  density  function,  p(n),  of  wave 
amplitude,  n»  of  the  form 


p (n)  = 


/2ifa 


exp 


HS)  ■ 


(4) 


for  different  currents  and  wind  speeds, 
b.  Wave  Frequency  Spectrum 

Hasselmann,  et  al.  (1973)  have  proposed  a one-dimensional 
fetch- limited  wave  frequency  spectrum  for  wind-generated 
waves  from  the  results  of  their  studies  in  JONSWAP  I which 
were  conducted  at  the  same  site  as  that  of  JONSWAP  II.  This 
wave  frequency  spectrum,  4>Q(n),  has  the  form  of 


4>o (n ) “ exP  t”1*25  (nm/n)4] 

n 


(5) 


where  nm  is  the  wave  frequency  at  the  peak  of  the  spectrum, 
K is  a parameter,  and  G(n)  is  a shape  function  which  is 
defined  as 


G(n)  - Yexe  ‘-(n  - V2/^2  %21 


(6) 


These  four  parameters , nm,  K,  y and  £ in  Equations  (5)  and 
(6)  are  determined  by  the  technique  of  simultaneously 
optimizing  all  four  parameters.  They  are  evaluated  as 


nm  = 21.99  (g/W)  (gF/W2)"0,33  , 


(7) 


K = 0.4775  (gF/W2)”0,22 


(8) 
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I • 


Y » 3.3 


9 


(9) 


0.07 

i 

0.09 


for  n in  (o,  n ) 
m 


for  n in  (n  . °°) 
m 


(10) 


.j 


where  F is  the  fetch  in  meters  and  W is  the  wind  speed,  in 

m/sec,  measured  at  10  meters  above  the  mean  sea  level.  The 

value  of  n^  given  by  Equation  (7)  is  very  close  to  the  value 

calculated  for  n„  by  Chen  (1972). 

m 

Let  8 be  the  equilibrium  coefficient  as  defined  by 
Phillips  (1958),  then,  from  Equations  (5),  (6),  (7),  (8)  and 
(9) 


0.4775 


\W  / 


-0.22 


3.3 


exp 


nW  / gF  \ 

' . * V ) 


0.33 

-21.99 


/967.22£‘ 


(11) 


One  can  notice,  from  Equation  (11) , that  8 is  not  only  a 
function  of  the  dimensionless  fetch  parameter  gF/W  but  also 
it  is  a function  of  the  dimensionless  wave  frequency  parameter 
nW/g.  For  fixed  values  of  W and  F,  as  the  value  of  n 
increases  the  value  of  8 approaches  a limit  which  can  be 
expressed  as 


8 = 0.4775  (gF/W2)  “°*22 
n <» 


(11) 
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A2 . NUMERICAL  RESULTS 

Theoretically,  wave  spectra  and  wave  height  distribu- 
tion functions  for  wave  interacting  with  various  steady  non- 
uniform  currents  can  be  evaluated  for  the  locations  at  which 
the  four  files  were  taken  on  September  19.  With  the  use  of 
Equations  (2)  and  (5),  wave  spectra  can  be  evaluated  by 
considering  wave  spectra  in  the  form  of  Equation  (5)  as  the 
wave  spectra  of  waves  without  the  influence  of  currents.  By 
the  same  token  wave  height  distribution  functions  can  be 
evaluated  by  using  Equations  (4)  and  (5).  Figures  A-l,  A-2, 

A- 3 and  A-4  show  the  theoretical  wave  spectra  of  waves 
interacting  with  currents  at  the  locations  where  File  4 of 
Tape  1,  and  Files  1,  2,  and  3 of  Tape  2 were  taken,  respec- 
tively. The  vertical  coordinates  are  presented  in  logarithmic 

2 

scale  and  in  the  unit  of  m /Hz  and  the  horizontal  coordinates 
are  presented,  also,  in  logarithmic  scale  and  in  the  unit  of 
Hz.  As  shown  on  the  top  of  these  figures,  the  currents' 
speeds  for  the  curves  from  the  top  down  are  -2,  -1.5,  -1, 

-0.5,  0,  0.5,  1,  2,  3,  and  5 m/sec  in  an  increasing  order. 
Figures  A-5,  A-6,  A-7,  and  A-8  are  the  theoretical  wave 
height  distribution  functions  of  waves  interacting  with 
currents  at  the  locations  where  File  4 of  Tape  1,  and  Files 
1,  2,  and  3 of  Tape  2 were  taken.  The  current  speeds  for 
the  curves  from  the  top  down  along  the  vertical  axes  are  -2 , 
-1.5,  -1,  -0.5,  0,  0.5,  1,  2,  3,  and  5 m/sec.  The  computer 
program  used  to  compute  for  Figures  A-l  to  A-8  is  included 
in  this  report  as  Appendix  B. 
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Theoretical  wave  spectra  of  waves  interacting  with  currents  at  the  location 
where  file  4 of  tape  1 was  taken 
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HAVE  FREQUENCY  HZ 

Fig.  A-2  - Theoretical  wave  spectra  of  waves  interacting  with  currents  at  the  location 

where  file  1 of  tape  2 was  taken 


MfiVE  FflEOUENCT  HZ 

Fig.  A -3  - Theoretical  wave  spectra  of  waves  interacting  with  currents  at  the  location 

where  file  2 of  tape  2 was  taken 
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^*8*  A-5  - Theoretical  wave  height  distributions  of  waves  interacting  with  currents 
at  the  location  where  file  4 of  tape  1 was  taken 
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Wave  Amplitude 

ftieoretical  wave  height  distributions  of  waves  interacting  with  currents 
at  the  location  where  file  1 of  tape  2 was  taken 
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MULTIPLES  QF  .810  METEHS 
Wave  Amplitude 

Fig.  A-7  - Theoretical  wave  height  distributions  of  waves  Interacting  with  currents 
at  the  location  where  file  2 of  tape  2 was  taken 
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